On the gravitational origin of the Pioneer Anomaly 
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From Doppler tracking data and data on circular motion ol astronomical objects we obtain a 
metric of the Pioneer Anomaly. The metric resolves the issue of manifest absence of anomaly 
acceleration in orbits of the outer planets and extra-Pluto objects of the Solar system. However, it 
turns out that the energy-momentum tensor of matter, which generates such a gravitational field 
in GR, violates energy dominance conditions. At the same time the equation of state derived from 
the energy-momentum tensor is that of dark energy with w — —1/3. So the model proposed must 
be carefully studied by "Grand-Fit" investigations. 
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I. INTRODUCTION 

Spacecrafts Pioneer 10 and 11 were launched in the 
early 1970's for the exploration of outer planets of the So- 
lar system (see the special issue of Science 183, No. 4122, 
25 January 1974, especially [H, 0]). After the encounters 
with Jupiter and Saturn they followed hyperbolic trajec- 
tories on leaving the Solar system. Because of rotational 
spin-stabilization of these spacecrafts, which reduces the 
need for manoeuvres, they represent unique experiments 
for testing our understanding of celestial mechanics. The 
accuracy of acceleration measurements for the Pioneer 
spacecrafts is about 10"^" m/s^ [1, 

During the flight spacecrafts were continuously tracked 
by Doppler effect on retransmitted radio signals. Then 
data were fitted to theoretical ones obtained in PPN ap- 
proximation initially by ODP program of JPL (JPL's 
Export Planetary Ephemeris DE405 was used for planet 
motion) . 

But surprisingly above 10 a.u. of heliocentric distance 
the systematical deviation of experimental and theoreti- 
cal data was found [5] . This deviation can be described 
simply as a constant acceleration towards the Sun with 
magnitude of about 8-10~^° m/s^. This value is the same 
— within error limits — for all the spaceships Pioneer 10 
and 11, Galileo and Ulysses and for all distances from 
the Sun 0. 

This coincidence has been interpreted as a hint of the 
gravitational — metric — origin of the acceleration. But 
at the same time there are no signatures of such an ac- 
celeration in the orbits of outer planets and other objects 
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in the Solar system. Inclusion of such acceleration leads 
to unavoidable deviations from the observed planet posi- 
tions 0,11. 

Many attempts to explain the anomaly were made dur- 
ing last 10 years. Some of the recent work on this subject 
includes analyses of: the thermal radiation of the Pio- 
neers l9l lOl. the gravitational attraction by the Kuiper 
Belt fill 12| . [l3ll . the cosmological origin of the Anomaly 
[3, llSl. Il6l [iTI. [Tsl. the influence of multipole mo- 
ments of the Sun [20|, the clocks accelerati on l2ll| and 
many proposes of modified gravity [H, [H, [13, |25| includ- 
ing even laboratory investigations on very small accel- 
eration dynamics [26| and interesting endeavors to con- 
strain some parameters of modified gravity theories by 
the known value of the Pioneer Anomaly [13, [1^ . 

In the frame of metric theories of gravitation there is 
an attractive possibility to explain the Pioneer Anomaly 
by metric perturbation, preserving at the same time the 
character of planet motion. It is possible because the 
Pioneers' trajectories are very different from planet or- 
bits: spacecrafts leave the Solar system almost in a radial 
direction, while planets orbit the Sun almost circularly. 
The potential possibility of such an approach was noted 
independently in a series of manuscripts [1^, [sO, [sTj , but 
authors of these works didn't analyze the origin of the 
metric. Our approach is closer to that of Kjell Tangen 
js^ . but with more rigor because we don't neglect per- 
turbation of the space components of metric. 

The goal of paper is to find the static space-time metric 
close to Schwarzschild one, in which radial motion of test 
bodies shows the Pioneer Anomaly, but circular motion 
doesn't. For this purpose in section[n]we develop and dis- 
cuss an algorithm of metric determination from data on 
radial and circular motions. Metric determination does 
not make use of Einstein equations so it is applicable to 



2 



any pure metric theory of gravity in terminology of Will 
et. al. [H, [131 • Then we apply the method in the case 
of the Pioneer Anomaly, starting from the Schwarzschild 
space-time and recover the properties of matter forming 
such a metric within GR (section HIH) . Finally some con- 
cluding remarks are made. 



II. SPACE-TIME DETERMINATION FROM 
RADIAL AND CIRCULAR MOTIONS 

A. Metric choice and time definition 

We begin with the interval with 3 metric functions r, 
p and a (for simplicity it was taken c = 1) 

ds^ = e^^^de - eP^'-^dr^ - e<^\\de^ + cos^ Od^^), 

(1) 

and then find the gauge relation between them for max- 
imal simplification. The metric functions r, p and a will 
be referred as time, radial and transverse metric func- 
tion or coefficient, respectively. We consider radial and 
circular motions in such a space-time separately and find 
connections between metric functions following from the 
known properties of the motion. But first of all wc must 
recall some time convention. 

The metric ([1]) is written in the form that is consis- 
tent with global clock synchronization, in fact d/dt be- 
ing timelike Killing vector. So in the approach proposed 
when the cosmological effects is totally neglected and So- 
lar system is supposed to be placed in a space-time, that 
is Minkowskian at spatial infinity, the coordinate time 
is the astronomical ephemerides time ET (up to a mul- 
tiplier). This is the usual approximation used for PPN- 
ephemeride calculations. 

There is a difficulty, because the perturbations re- 
quired by the Pioneer Anomaly grow with the radial dis- 
tance, so the perturbed space-time is not asymptotically 
flat. But it is not a big problem because before the met- 
ric perturbations grow significantly the space-time has 
a wide nearly-flat region in which we can choose almost 
Minkowskian observers and coordinates. So further if we 
talk about " spatial infinity" , we mean this wide region, 
in the Solar system scale the effects of difference between 
this approximation and rigorous treatment are negligible. 
The same problem with the same solution is arising in 
the PPN-approximation then we must place the system 
not in Minkowskian background but in the cosmological 
one. Additionally, as it can be shown, in PPN-approach 
the cosmological effects, such as mutual acceleration of 
geodetically moving bodies, have the second order in H . 
Therefore even while the Pioneer Anomaly acceleration 
is nearly equal to ci?, in the framework of pure metric 
theories of gravitation there is no possibility to link it 
to the cosmological expansion (early but almost exhaus- 
tive analysis of the problem was made by R. C. Tolman 
[ssl . §§153-156], for the recent work on subject see the 



articles mentioned in the Introduction and additionally 
[H, [lljls^ and references within). 

The radial coordinate rescaling r — > /(r) changes all 
metric functions, giving a possibility to imply gauge con- 
ditions on the metric coefficients. But there are two 
invariants, i.e. physically measurable quantities, which 
characterize the distance from a given point to the cen- 
ter of space symmetry. Firstly it is an atomic time rate in 
comparison with the coordinate time rate (or atomic time 
rate on "spatial infinity") e"^/^, and secondly it is an area 
of a sphere of points equidistant to the center of the space 
47rr^e'^. So the numerical values of time and transverse 
metric coefficients have clear physical meaning for a given 
space-time point and only p{r) can take arbitrary values. 
Usual choices include p = a corresponding to isotropic 
coordinates of PPN-approximation and implicit on p(r) 
relation ct = corresponding to Schwarzschild coordi- 
nates. 



B. Radial motion and its description by Doppler 
radio tracking 

Radial motion in the space-time is fully described by 
the energy gttjf = e'^*-''-'u° = k = const and the 4- 

velocity length conservation e'^^^'-'u"^ —e''*^'"'?/^^ = e, e = 
or 1 for electromagnetic waves and test bodies, respec- 
tively: 



dr ^l-ee^W/fc2' 

the constant k being connected with the velocity v of the 
spacecraft on "space infinity" 

= (3) 

The relation ^ can be integrated to give us t{r) depen- 
dence which in turn can be inverted giving r{t). But 
there is a gauge freedom in the result because we can 
choose p{r) freely. Moreover, this relation essentially in- 
volve both physical r and unphysical p. 

Our goal is to determine the metric functions from ob- 
servations. Direct results of experimentation in astron- 
omy and cosmology are the measurements of externally 
originated signals received on the world line of an ob- 
server. These signals have mostly electromagnetic char- 
acter and can originate from some external source (e. g. 
the light emitted by the Sun and reflected by a planet) or 
be emitted by the observer him/herself and then returned 
to him/her after some interaction with outer bodies (as 
in case of radar measurements). The latter case is con- 
sidered here as corresponding to the real situation jSj . 

The scheme of Doppler tracking used in the Pioneer 
experiments is very simple conceptually: an electromag- 
netic signal, emitted from the world line of the observer, 
is reflected back by the mirror on the world line of the 
spacecraft and then compared with the initial one again 
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on the observer world line. To be more precise, the 
monochromatic electromagnetic signal, obtained from 
the high precision hydrogen maser 0, subsection III. A], 
is emitted by the antennae on Earth in the direction of 
the spacecraft. This signal is detected by the spacecraft, 
amplified and reemitted back to Earth, where the mea- 
sured waveform of arrived signal is compared with the 
emitted one to obtain the red shift of the signal and the 
time of signal travel (the very detailed description of the 
process can be found in Q). 

Now we must describe the Doppler tracking and the 
signal time arrival analysis in this space-time. In the 
geometric optics approximation (which is applicable for 
the case considered) the Doppler shift is governed simply 
by the ratio between scalar products of the 4- velocity on 
world lines and the null wave vector of the signal, parallel 
transported along the null geodesic line between emitter 
and receiver: 



u„ ■ kp 



(4) 



where Vj. and v^, are received and emitted frequencies, 
measured by the standard atomic clocks, 
Sr and Se are proper times of one cycle of oscillation, 
Ur and Ue are 4-velocities of receiver and emitter, 
kr and ke are tangential null vectors (wave vector), par- 
allel transported along the path of the signal. 

Atomic clock time deviations from ephemerides time 
along with all known effects of Earth motion were taken 
into account during the data processing (see Q), so for 
the description of such a small deviation like the Pioneer 
Anomaly it is sufficient to use the simple model, in which 
the emitter of the initial signal and the receiver of the 
retranslated signal are fixed at constant distances from 
the Sun on the line from the Sun to apparatus. The 
signal is emitted from this "fixed" Earth at rg and t — tp, 
received by the spaceship at r and t, amplified, exactly 
retransmitted back to Earth and finally compared with 
the initial frequency on the " fixed" Earth again at rg and 
t + tp {tp is the time of signal propagation, the same for 
forward and backward directions). 

As it can be shown easily, the frequency Vr received on 
Earth is connected to the initially emitted v^ as 



1 - Vl - e^W/fc2 
''1 + Vl - e^W/fc2' 



(5) 



This expression can be readily reduced to special rela- 
tivistic one in the case of e'^ = 1. 

As a relation between physical quantities only, this 
equation does not involve arbitrary unphysical p func- 
tion. The problem of gauge choice comes with the def- 
inition of r{t): from ([5]) one can determine e'^^'''-*-*\ but 
because the radial coordinate r (and consequently r{t)) 
is arbitrary, the radial dependence of the time metric co- 
efficient remains gauge dependent. It is interesting also 
that this relation cannot be represented by a power series 
in terms of small deviations of e'^^-^^ and k from 1. The 



transverse space metric coefficient a(r) naturally cannot 
be determined from the radial motion only. 

So we come to an unavoidable alternative of a-priory 
r(t) definition or a-priory imposing some gauge condition 
on p{r). In each case the remaining function is defined 
by the experimental data, and our goal now is to find in 
which case the process of metric restoration can be done 
without unnecessary complications. In the next subsec- 
tion we consider these possibilities in some details and 
then show that the best choice is the latter case, i.e. im- 
posing a gauge. 



C. General formulae and coordinate choice 

While the time of signal arrival to the spaceship t can 
be easily determined as a half-sum of observed times of 
emitting and receiving tr of the signal at the "fixed" 
Earth 



t 



(6) 



the corresponding r determination is not so trivial task. 
In general we can measure only the time of signal travel 
from the "fixed" Earth to the spaceship as a half- 
difference between observed times of sending and receiv- 
ing of the signal 



tr 



lM tr — tf, 

ar = — - — 



(7) 



These are results of a different method of tracking — 
signal time arrival analysis, which in essence represents 
integration of the Doppler data. 

So to recover p{r) from a given r(t) we must ffistly find 
T(r) by ^ from the observed redshift, and then solve an 
integral equation above. On the other hand, to find r(t) 
from a given p{r) it is sufficient to solve a non-integral 
equation following from relation ^ for the spacecraft 
motion 



to 



e 2 ar, 



(8) 



where to is the time when the spacecraft leaves "fixed" 
Earth. The left-hand side of the relation can be found 
totally from the observed frequency shifts by ^ , and the 
right-hand side is a known function of r with given p{r). 

Consequently maximal simplification of the problem is 
reached in the case of a-priory given radial metric func- 
tion p{r). It is nonsense to define it dependent on still 
unknown time and transverse metric functions, with one 
interesting exclusion: if r(r) = p{r) then r(t) can be 
recovered from ([7]) simply as 



(9) 



This choice of coordinates known as light or null coordi- 
nates is not so usual as Schwarzschild (tr = 0) or isotropic 
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{p = a) coordinates but it is the most suitable one for the 
situation. Both abovementioned choices are especiaUy in- 
appropriate here because they rely on transversal metric 
function that does not reveal itself in pure radial motions. 



D. Circular motion 

We know that the near-circular motion of outer Solar 
system objects (i.e. Neptune or PlutoHs unperturbed by 
the Pioneer Anomaly acceleration This gives us a 

way to determine the transversal metric coefficient and 
therefore to find metric completely. 

The angular velocity of circular motion (0 = 0, </> = uot) 
in the considered space-time is defined by the ratio of 
derivatives of time and transversal space metric coeffi- 
cients 



(10) 



where as usual prime denotes differentiation with respect 
to radial coordinate, and dot denotes differentiation in 
time. So there is no dependence on the radial metric 
part at all. Moreover as uj {r{t)) and (e^('^(*»)' is directly 
observable, so the transverse metric coefficient r^e'^'-''*^*-'-' 
can be obtained by a simple integration without any no- 
tion of the radial metric function. 

So again recalling the radial motion equation ([2]) we 
conclude that all gauge conditions for p{r) involving 
transverse metric coefficient are not convenient for treat- 
ment of the Pioneer Anomaly, because all such conditions 
lead to coupling of equations ([5]) and (|10p which can be 
solved independently otherwise. 



E. Final list of relations and concluding remarks 

So we work in light or null coordinates T(r) = p(r), 
then: 



dt 



1 



dr ^1 -ee^W/fc2 



(11) 



The trajectory of the spacecraft r{t) is recovered simply 
from time of sending te and arrival tr of signal 



in 



r - To 



t 



ijf I i g 



r = ro + 



and the observed redshift of the signal 
z{t) = 



I/e (l-e^('-(*))/fc2)-l/2 + l 



(12) 



, (13) 



can be immediately transformed into time metric coeffi- 
cient 



oT(r(t)) 



2 - z{t) 



{2-zy 



(14) 



The transverse space metric coefficient is defined by the 
dependence of angular velocity uj on radial coordinate r 



Ak'^z{r)z'{r) 
(2-z(r))3tj2(^) - ' 



dr, 



(15) 



which for the space-time of the Pioneer Anomaly should 
be the same as in the Schwarzshild field (or very close 
to it). So basing on the Schwarzshild metric one can 
find such perturbations of metric coefficients, that the 
circular motion remains unperturbed, but the radial one 
shows small deviation — the Pioneer Anomaly. 

It is worth noting that the sentence "the circular mo- 
tion remains unperturbed" denotes exactly the following: 
if, neglecting all mutual planet disturbances, the period 
of circularly orbiting planet will be measured in the units 
of ephemerides time and the distance of planet from the 
baricenter of the Solar system will be found by analyzing 
light propagation times on straight lines from the Sun 
(null coordinates!) then the values of the period and the 
distance will exactly be the same as needed for the 3rd 
Kepler law to hold — exactly as in the Schwarzshild field. 



III. THE PIONEER ANOMALY AND ITS 
SOURCE IN GR 

A. Schwarzschild space-time in light coordinates 

Radius r in light coordinates of the Schwarzschild field 
is related to the usual Schwarzschild radial coordinate 

as 

r = + rg In f ^ - 1 j , rs = Vg {l + W [e"^ 



(16) 



The interval in null coordinates is as follows 



(17) 



so that 



X^=^^.W/;(e--M, (18) 



l + W 



-W = J(l-,M.(e^-))^ 



(19) 



where W{x) is the so called multiplicative logarithm or 
Lambert W function 



T^(a;)e^(") = X. 



(20) 



B. The Pioneer Anomaly. Radial perturbation 

The Pioneer Anomaly is the hnear in ephemerides time 
ET deviation of the experimentally obtained frequency of 
received signal from the "modelled" one 



2a p 



(21) 



where ap 
eration 16 



8 • 10^^" ni/s^ is the "unmodelled" accel- 
The ephemerides time coincides with time 
coordinate t of the metric considered (as described ear- 
lier). For the most part of range of the Pioneer Anomaly 
found (from ~ 15 to ~ 80 a.u.) the deviation of the Pio- 
neers' orbits from pure radial motion is comparable to or 
even below experimental uncertainty in the acceleration: 
it can be estimated roughly as a ratio of the semi-major 
axis a absolute value to heliocentric distance r 



< 



lO^km 



17%, 



(22) 



r ~ 40 • 150 • lO^km 
while experimental error in the acceleration is 
1.33/8.74 ~ 15% (see Appendix of 0). It should 
be noted that this uncertainty prevents Anderson et. al. 
from determining the direction of the acceleration: to 
Earth or to the Sun (see beginning of section VII and 
especially note 73 of @). 

The modelled frequency and velocity of the spacecraft 



1 - 



\ 



1 - 



1 -i;2 



|) =...„(.■) = 



\ 



l/W{e 



(23) 



(24) 



so accurately expanding the expression (fH|) for the time 
metric coefficient from the red shift z{r{t)), one can find 
that given the accuracy of the Pioneer Anomaly measure- 
ments one can simply use the relation 



z{r) Sz{r) f \A+f \ 



(25) 



where At = t{r) — to is the time from the start of the 
Pioneer Anomaly (we suppose that for r < Tq the metric 
coincides with the Schwarzschild one, so before to = t(ro) 
there is no anomalous acceleration), Sz{r) is a deviation 
of observed red shift from the modelled one. 

In the first approximation t{r) dependence can be re- 
placed with the "modelled" time, which to the accuracy 
of the measurements is the same as in the Newtonian 
case 



tniir) = t{ro) 



ro 



r 



— ^ sinh 



—V . (26) 




r, a. u. 



70 



FIG. 1: Metric perturbation Se'^^^^ of tlie Pioneer Anomaly for 
V from 5 km/s to 50 km/s in 5 km/s steps (solid lines from 
bottom to top) for the metric matching the Schwarzschild 
metric at 12 a. u., compared to the "naive" post-Newtonian 
one (dashed line) 



Finally inserting this and modelled z{r) into the equation 
(|25p we arrive to the perturbation of the time metric 
coefficient 



(5e^('') = 2af 



-C 



1- Wl 



ZE^sinh-M./^z; 



(27) 



where C is a constant which can be determined from ro 
and V. The perturbation appears to be non-linear in r, 
but for the Pioneer 10/11 parameters the deviation from 
linearity is buried deep in the experimental errors. It is 
illustrated by the figure [l] which shows the deviation of 
the time metric coefficient compared to the "naive" post- 
Newtonian approach, where one simply adds to the grav- 
itational potential <I>(r) a term linear in radius and use 
^■r{r) ^ I ^ <I>(r). As we can see, the difference is mainly 
in the slope, all the graphs are nearly linear. Moreover, 
the relative value of the deviation from linearity is de- 
creasing with radial distance. So the linear approxima- 
tion Je"^^''-' ~ 2riap{r — ro) is sufficient for the Pioneer 
Anomaly explanation. The only difference between this 
more accurate result and the " naive" post-Newtonian one 
is the presence of rj, which is always less then 1 (see table 



C. Transversal perturbation leaving planet orbits 
unchanged. Necessity of the central source 

Now we assume that in the perturbed space-time plan- 
ets orbit with the same periods as in the unperturbed 
Schwarzschild solution, so that there is no signatures of 
the Pioneer acceleration in the orbits of planets. So the 
angular velocity dependence uj{r) must be the same as in 



6 



TABLE L The quantity 1 — r; of the best hnear approxunations ^e^''^' ~ 2r)ap(r ~ ro) on the interval ro < r < 70 a. u. for 
different velocities v and metric matching distances ro* 



ro, a. u. 










V 


, km/s 










5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


10 


0.175 


0.098 


0.057 


0.036 


0.025 


0.018 


0.013 


0.010 


0.008 


0.007 


15 


0.146 


0.077 


0.044 


0.027 


0.018 


0.013 


0.010 


0.008 


0.006 


0.005 


20 


0.122 


0.062 


0.034 


0.021 


0.014 


0.010 


0.007 


0.006 


0.005 


0.004 



* ro is the radial coordinate at which the metric coefficients coincides with the 
Schwarzschild ones. 



the Schwarzschild case (see, e. g., eq. (25.40) of [3 



1 



2rl{l + W {e^o ^ 



2rV 



(28) 

and the perturbation in must lead to such a pertur- 
bation in r^e"^ that uP'ir) remains invariant. In the per- 
turbed space-time by the general formula pop we have 
from simple mathematics 

2n = ^I^= (e-('-) + fc-('-))' 



(e^W)' + (5e^W)' 



(5e"W)' 



(r2e'^M)' + ((5(r2e'^('')))' ((5(r2e'^W))' ' ^^^^ 
So the perturbation of the transverse metric coefficient is 

(^(r^e'^M) = r^Je'^W = /" Lo-'^{r)5{e''^'''^)' dr, (30) 

J ro 

and finally 

3 



r 

iapr) r 



ro 

2 „r / \ 3 
5 



r 



ro \' 9 



dr 



dr ~ 
4apri{r'^ - r^) 



(31) 



In the first approximation the perturbation of r^e"''^''^ 
grows quartically in radius. 

Note the gravitational radius of the source rg in the an- 
swer. So the effect of the Pioneer Anomaly can he repro- 
duced only by perturbations of Schwarzschild space-time 
and not Minkowski one. So the gravitational explana- 
tion of the Pioneer Anomaly can be obtained without 
the equivalence principle violation required by various 
authors 0,1,112. 



D. Matter corresponding to the obtained metric in 
General Relativity 



lution found for the weak field of a point mass. But we 
think that it is not very promising because there are no 
reliable experimental evidence in favor of any gravita- 
tional theory other than General Relativity. 

Instead we find the properties of matter surrounding 
the point mass in GR that can generate the obtained 
metric. Because in the scale of Solar system experiments 
the influence of cosmological constant is negligible, for 
the determination of matter one can use the Einstein 
equations of the form 



Gij — R 



SttG 



(32) 



Using for simplicity the metric in the form 

ds^ = e^^^'^dt^ - eP^'''Ur'^ - e'^^^'^dO^ + sin^ 9d(j)^), (33) 
one arrives at the Einstein tensor 



^ — i^XtT, ® T, - A,5, Sj - j , (34) 

5, - {0,e^0,0}, = {e^,0,0,0}, (35) 

-S^S' = T,T = 1, S'T, = 0, (36) 

At = AeP-" + {p' - 2a' - r') (a' -t') + 2 (t" - a") , 

(37) 

A, ^ AeP--" - t' {a' - r') - p' {a' + t') + 2 (t" + a") , 

(38) 

A = + cj't' + r'2 - p' {a' + r') + 2 (r" + a") . (39) 
Expanding to the first power of ap one obtains 



At = -96 , \s 

^9 



-32^4. 
rg r4 



(40) 



The algebraic type of the energy-momentum tensor at 
spatial infinity is that of an ideal fluid (by A^ 0) with 
constant positive pressure 



One can try to find the gravitational field theory, that 
gives equations of the gravitational field allowing the so- 



e , apr\ 
p^ — A-> 12— > 0, 



4k 



nra 



(41) 
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but negative energy density 

p^^—[\,-X)^-m^<Q. (42) 

4 k KVg 

It is worth noting that relation between p and p is as 
for an uhrarelativistic fluid except for the sign: instead 
of p — p/3 one has asymptoticaUy p = — p/3. This is a 
typical equation of state of dark energy with parameter 
w = —1/3. It is interesting that such a fluid does not 
change the cosmological dynamics of Friedman-Lemaitre 
universe (see, e.g., 140|, III.E]). 

IV. CONCLUSIONS 

In this work we show that it is possible to perturb time 
metric coefficient of the Schwarzschild space-time in such 
a way that the Pioneer Anomaly is reproduced. More- 
over, because planet motion is governed also by anther 
component of the metric, we can tune it so that circu- 
lar orbits is not disturbed. This result applies in any 
pure metric theory of gravitation where test bodies follow 
geodetics of the metric. It is deduced that the perturba- 
tion of the time metric coefficient can be taken linear in r 
without contradiction with the accuracy of experimental 
data obtained up to now. 

Assuming the validity of the General Relativity, we 
find out the energy-momentum tensor generating the 
metric obtained. The tensor corresponds to an ideal fluid, 
however having negative energy density. It is interesting 
to note that an exact static solution with spherical sym- 
metry is known for the "fluid" with ip-\- p = Q [4l|, §8.5]. 
This "fluid" does not interact with the ordinary matter 
besides its gravitational influence on the metric, so it 
much like WIMPs or scalar field of gravitational theories 
of Brans-Dicke type. Thus the ordinary matter including 
spacecrafts and planets is moving geodesically. 



Naturally the found "fluid" does change the planet or- 
bits (if they are not strictly circular) and light rays paths. 
The model proposed must be carefully studied in view of 
the " Grand-Fit" investigations [i^, , but direct mea- 
surements from the planned missions for testing General 
Relativity in space are preferable [3, E, . The ab- 
solute value of effects for the perturbation found as well 
as for exact solution of Stanjukovich and Ivanov will be 
studied in the forthcoming paper. 

The analysis presented in this paper can encourage 
someone to find out which of the known alternatives to 
GR can reproduce the metric found or to invent some new 
theory which can do it. It is possible, but in our opin- 
ion the Pioneer Anomaly has some simple explanation by 
conventional and non-gravitational physics, which is not 
found yet. So the question of deep theoretical grounds 
for the existence of the "fluid" in GR or of development 
of some new theory of gravitation based on the metric 
obtained is not in the scope of our article. Instead we 
point out that negative energy density of the "fluid" is in 
a direct contradiction with the properties of conventional 
matter. One interesting possibility for the "fluid" is dark 
energy, which has the right equation of state. Therefore 
we suppose that at present the metric (gravitational) ori- 
gin of the Pioneer Anomaly cannot be ruled out. 
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